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Abstract 

The effects of the ring-exchange Hamiltonian H3 — J3 ^^(^^^^(S'i • Sj){Si ■ Sk) on the triangular lattice are studied using classical 
Monte Carlo simulations. Each spin Si is treated as a classical XY spin taking on Q equally spaced angles (Q-states clock model) . 
For Q = 6, a first-order transition into a stripe-ordered phase preempts the macroscopic classical degeneracy. For Q > 6, a finite 
window of critical phase exists, intervening between the low-temperature stripe phase and the high-temperature paramagnetic 
phase. 
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The importance of higher-order ring-exchange processes 
in low-dimensional magnets and its potential role in stabi- 
lizing liquid-like exotic ground states has recently been un- 
der intense investigation[l,2]. Generically the Hamiltonian 
is of the type 



(1) 



Tl>2 



^ ■ with Hn for n ~ 2, 3, 4 given by 



H2 — J2 Si ■ Sj 

{i]k) 

Hi — J4 ^ {Si ■ Sj){Sk ■ Si). 

{ijki) 



(2) 



Each Si is a Heisenberg spin, and {ij), (ijk) and (ijkl) 
refer to nearest-neighbor pair, triplet, and quartet of sites, 
respectively. 

In particular the possibility to realize a stable spin-liquid 
ground state due to H3 in a two-dimensional triangular 
lattice has been suggested in the variational Monte Carlo 
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study of Motrunich[2]. It may be inferred that the three- 
site exchange process has the "frustrating" effect which 
renders the liquid ground state energetically more stable 
over a magnetically ordered structure. 

In this paper, we report the first results of isolating the 
effects of the three-site exchange process by studying the 
following model. First, we consider the classical counter- 
part of the Hamiltonian (1) where Si is treated as a uni- 
modular vector. Secondly, only H = H3 is considered in 
this paper while leaving the study of the compound models 
such as H = H2-\- or H = H2 + i^a + for the future. 
Thirdly, we consider the planar spin Si — (cos^i, sin^i). 
The angle 9i is divided up into Q equally spaced segments. 
The same strategy had been applied for H = H2 as a way 
to asymptotically approach the behavior of the XY model 
{Q ~ 00) and is known as the Q-states clock model[3,4]. 

Antiferromagnetic Q-states models on a triangular lat- 
tice have double transitions of XY- and Ising-types with 
extremely close critical temperatures [5]. The same sub- 
tlety might pervade the H = H3 model too, but here we 
choose to focus on the broader issue: What is the nature 
of the low-temperature phase exhibited by H = H3 ? The 
results for Q = 2 and Q — 6 are discussed in detail in this 
paper. Some preliminary specific heat data for larger Q 
are presented. 

Q=2: It turns out that Q = 2 model maps onto the 
antiferromagnetic Ising model on the triangular net, first 
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studied by Wanmer[8]. In the Ising case spins take on Si 
±1, and the Hamiltonian reduces to 



{Si ■ Sj){Si ■ Sk) SjSk, -ffs "^Ji 



3 ^ SiSj, 



(3) 



which is the antiferromagnetic Ising model. This model pos- 
sesses macroscopic degeneracy[8,9] which is also revealed 
as the residual entropy Sq. Our Monte Carlo (MC) calcu- 
lation gives So ~ 0.323fcB, in excellent agreement with the 
value predicted earlier[8,9]. There is no long-range order 
down to zero temperature in this model. 

With Q > 3 H2 and H3 are no longer equivalent. The 
lowest-energy configurations for consist of two-up and 
one-down (or vice versa) spins for each elementary trian- 
gle. Thus, macroscopic degeneracy is a general feature of 
H = H3 for an arbitrary even integer Q. It is also well 
known that the magnetic ordering for H = H2 is obtained 
for angles of 120° between nearest-neighbor spins. Such 
situations are possible if Q is a multiple of 3. To allow the 
realization of both, we consider the case Q = 6. 

Q=6: With Q — 6 we observed a first-order transition 
at Tc/Js « 1.05. Hysteresis in the average energy and the 
order parameter (defined below) in our MC runs vindicate 
the first-order nature, as shown in Fig. 1. 
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Fig. 1. Low-temperature configuration for the Q = 6 model, H = H^. 
Ferromagnetic ordering within a diagonal stripe and antiferromag- 
netic ordering between adjacent stripes are apparent in (a). Both the 
average energy (b) and the average magnetization (c) data are con- 
sistent with the first-order phase transition to the low-temperature 
phase. 



The nature of the low-temperature, ordered phase is 
clearly demonstrated in Fig. 1. We find the spontaneous 
emergence of stripe-like domains of up and down spins be- 



low Tc- The ground state degeneracy is lifted through the 
order-from-disorder mechanism [6, 7]. While a conventional 
order-from-disorder idea predicts the gradual separation of 
the free energies of different classical ground state config- 
urations with rising temperatures and no phase transition, 
our model exhibits a first-order transition which pre-empts 
the classical macroscopic degeneracy. The difference is due 
to the discrete nature of our model. The order parameter 
appropriate for this stripe-like configuration is 



1 



(4) 



where each lattice site is given the coordinate i = 
iiBi + 1262, ei = X, 62 = —x/2 + \/3y/2, and N is the 
number of sites. 

Q>6: For a finer spin segmentation we still obtain the 
low-temperature stripe-like phase. A single first-order 
phase transition observed for Q = 6 is split into two tran- 
sitions, at temperatures Ti and T2 with T < Ti being the 
stripe-ordered phase. The intermediate phase Ti <T < T2 
appears to be critical[4,10]. 



- 48*48 

- 72*72 



Fig. 2. Specific heat C{T)/T for Q = 12, H = Hs, for 48 X 48 and 
72 X 72 lattices. The two peaks are indicative of the presence of two 
phase transitions. The low temperature phase is given by the stripe 
configuration shown in Fig. 1. 
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